Effects of spin fluctuations in the t-J model 
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Recent experiments on the Fermi surface and the electronic structure of the cuprate-supercondutors 
showed the importance of short range antiferromagnetic correlations for the physics in these systems. 
Theoretically, features like shadow bands were predicted and calculated mainly for the Hubbard model. 
In our approach we calculate an approximate selfenergy of the t-J model. Solving the U = oo Hubbard 
model in the Dynamical Mean Field Theory (DMFT) yields a selfenergy that contains most of the local 
correlations as a starting point. Effects of the nearest neighbor spin interaction J are then included in 
a heuristical manner. Formally like in J-perturbation theory all ring diagrams, with the single bubble 
assumed to be purely local, are summed to get a correction to the DMFT-self engergy This procedure 
causes new bands and can furnish strong deformation of quasiparticle bands. 



The obvious importance of antiferromagnetic 
(AF) correlations for the physics of cupratc su- 
perconductors found an impressing confirmation 
through ARUPS-experiments performed by Aebi 
et al. [§. There, in a Fermi surface (FS) map- 
ping of BiaSraCaCuzOslOOl) shadows of the FS 
were observed. The interpretation as an effect of 
short range antiferromagnetic correlations led to 
an intensive discussion of so called shadow fea- 
tures, which have been predicted by Kampf and 
Schrieffer § in a semi-phenomenological theory 
for the Hubbard model. 

In the paper we present an approximate eval- 
uation of the self energy for the t-J model Htj 
consisting of two parts: The hopping term 



Ht 



<1J>CT 



(1) 



acting only in the space of no double occupancy 
(c|g. = (1 — ni^_cr)c|^) is equivalent to the U — oo 
Hubbard model. The spin interaction is described 
explicitely by 



(2) 



<ij> 



The main idea of our approach is to take into 
account both, the spin fluctuations due to Hj 
and the strong local correlations in Hf. For the 
latter the dynamical mean field theory (DMFT) 
1^, ^, H , which becomes exact in the limit of infi- 
nite spatial dimensions, proved to be quite a good 



tool. It yields a strictly local selfenergy, which re- 
flects the metal-insulator transition for increasing 
U as well as a Kondo-like resonance around the 
Fermi energy. Even as an approximation for a 
two-dimensional system the DMFT-results seem 
to be reliable when being compared to e.g. Quan- 
tum Monte Carlo calculations ||, @, H) . We per- 
form the DMFT using the so called Non Crossing 
Approximation (NCA) ^, |l^ to solve the related 
local problem. For details of this procedure see 
ref. [|. 

The quantities obtained with this method 
are first the local selfenergy for calculating the 
approximate one-paricle Green's function of the 
twodimensional [/ = oo Hubbard model 



Gt(k, z) 
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and second the full local magnetic susceptibility 
X/oc('^)- The dispersion is given by the Fourier 
transform of the next neighbor hopping = 
COS kxd -\- cos kyo). 
In addition to the correlation effects already 
contained in Gt(k, z), spin fluctuations are now 
included by a formal perturbation expansion in 
Hj, taking Ht as the unperturbed part of the 
Hamiltonian and Gt(k, z) as the unperturbed GF 
(J = 0). 

As a first step all ring diagrams are summed 
resulting in a RPA-like structure for the effective 
interaction. In contrast to the standard RPA, 
however, the intermediate polarization diagram 



is the full susceptibility of the U = co Hubbard 
model. Since the latter depends only very weakly 
on q |pl| we further replace the latter quantity by 
Xiociz), which can be easily calculated from the 
effective local problem. 

This procedure results in an effective interac- 
tion 

V(q) = 3 J(q) 1 J(q) 

2 1 - J(q)xioc(ii^m) 2 1 + J(q)x/oc(it^m) ' 

from which we obtain a selfenergy contribution 
^fluctik) = -J2Gtik-q)V{q). (5) 

Here we use the notation k = (k, ia;„) and an 
implicite sum over bosonic Matsubara frequencies 
is contained in equations (4) and (5). E/;uct is 
calculated on the real frequency axis. The q-sum 
in eq. (5) is performed via Fast Fourier Transform 
on a 32x32 lattice and the sum over Matsubara 
frequencies by contour integration. 

The complete selfenergy is then given by 

i;(k, z) = ^lociz) + ^fluctO^, z) + ^Hartree, (6) 

where TiHartree IS a real number and here deter- 
mined by the occupation number of the DMFT- 
result. This scheme is iterated with the resulting 
Green's function replacing Gt in equation (5) un- 
til selfconsistency was reached. 

Fig. 1 shows the resulting spectral functions 
A(]s.,uj) (full lines) at the F— and the M— point 
(upper and lower picture), respectively for dop- 
ing i5 = 0.15, coupling constant J = 0.144t* and 
inverse temperature /3 = lO.Ot*. These two k- 
points would be equivalent in the case of antifer- 
romagnetic order. 

In comparison with the J = case (dashed 
lines) an additional peak above the Fermi energy 
in A(kr,w) is clearly visible, which we interpret 
as a shadow of the quasiparticle peak at the M- 
point. 

This becomes more clear in Fig. 2 , where 
the corresponding bandstructurc is plotted along 
high symmetry directions of the Brillouin zone. 
Crosses represent the maxima of A(k, u) when 
spin fluctuations are taken into account. The 
triangles correspond to the J — case. A new 
branch appears in the vicinity of the F-point 





0.45 














// 


\\ 
\\ 










// 


\\ 








0.30 


// 


\\ 










// 


\\ 










// 


\\ 






^< 




1 


\\ 








0.15 


// 












// 

// 




~ \ 






0.00 









-3 -2 -1 1 



3 



0.45 






!l ',1 
ll 


0.30 






1 ' 
I 1 " 

I H 


0.1 5 




/ 


1 \ 
\\ 






Jl 


\ 


0.00 







-3 -2 -1 1 

Figure 1: Spectral functions of the system with 
spin fluctuations (solid) and without (dashed) at 
the F- and the M-point. The solid upper curve 
(F-point) shows an additional peak, exactly at the 
position of the quasiparticle peak at the M-point. 
The parameters are doping S — 0.15, coupling 
constant J = 0.144t*, inverse temperature f3 — 

lo.or. 



above the Fermi energy, which shows the symme- 
try of the magnetic Brillouin zone, as is expected 
for shadow bands. As can be seen in both flgures, 
the main features due to strong local correlations 
are not altered by the spin fluctuation effects. 

In conclusion an extension of the dynamical 
mean field theory was presented including nonlo- 
cal spin-fluctuations in the t-J model. With this 
approximation we are able to obtain shadow fea- 
tures in the spectral functions. So far our numer- 
ical calculations yield shadow bands only above 
the Fermi surface. In the experiment ^ shadow 
bands cross the Fermi surface. We hope to repro- 
duce this efect by including a next-nearest neigh- 
bor hopping t' in our tight-binding bandstructurc. 

This work was supported by grant number 
Pr 298/3-1 from the Deutsche Forschungsgemein- 
schaft. 
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Figure 2: Bandstructure for the system with 
(crosses) and without (triangles) spin fluctua- 
tions. Around the F-point a new branch is visible, 
which shows the symmetry of a reduced Brillouin 
zone, although the system is paramagnetic. The 
parameters are the same as in Fig. 1. 
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